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!.  INTRODUCTION 


i he  Ballistic  Research  Laboratory  (BRL)  of  the  U.S.  Army  Is  designing 
a  Large  Blast/Thermal  Simulator  (LB/TS),  which  Is  a  proposed  facility  for 
creating  pressure  and  thermal  transients  like  those  that  occur  near 
nuclear  explosions.  (l} '  The  LB/TS  Is  basically  a  large  shock  tube  with 
area  changes  along  its  length.  Unlike  a  conventional  shock  tube,  which 
uses  rupture  diaphragms  to  separate  the  high  pressure  driver  region  from 
the  low  pressure  driven  region,  It  Is  desirable  to  utilize  fast-opening 
valves  in  the  LB/TS  for  wave  shaping  purposes  as  well  as  operational 
efficiency.  However,  with  real  valves,  a  finite  opening  time  results  in  a 
compression  wave  with  a  finite  rise  time.  The  question  then  arises  as  to 
whether  the  compression  wave  will  steepen  to  an  adequate  shock  profile 
within  the  specified  wave  propagation  distance,  to  simulate  the  desired 
blast  wave.  In  support  of  this  design  effort,  the  Mathematical  Analysis 
Unit  at  the  Idaho  National  Engineering  Laboratory  (INEL)  has  done 
calculations  to  determine  the  maximum  throat  valve  opening  times  allowable 
for  the  creation  of  certain  shock  waves  In  the  LB/TS  test  section. 
Calculations  were  done  for  9  specific  cases  specified  by  the  BRL, 
representing  9  different  initial  conditions  In  the  shock  tube  driver.  All 
the  calculations  were  done  with  the  BRL-Q1D  computer  code^.  This 
report  summarizes  the  results  of  the  valve  opening  calculations. 
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2.  CALCULATIONS 


For  the  present  studies,  the  LB/TS  was  repress  ' 

shown  In  Figure  1  and  described  In  Table  1.  We  are  interested  in  tne 
Initial  magnitude  and  rise  time  (not  the  later  decay)  of  the  shock  at  tu 
test  section  103.5  meters  from  the  valve.  The  objective  of  the  present 
study  Is  to  determine  how  fast  the  valve  must  open  w  a1 low  the 
overpressure  to  rise  from  5%  to  95%  of  the  plateau  value  in  0.5 
milliseconds  or  less.  We  are  not  concerned  with  reflections  from  tne  e~d 
of  the  shock  tube,  which  arrive  too  late  to  affect  the  shork  front  of 
Interest.  This  allows  us  to  rjdel  a  shock  tube  substantially  shorter  t;  :  > 
the  real  LB/TS,  thereby  saving  some  computer  costs. 

Calculations  were  done  for  the  9  cases  listed  in  7 able  2,  represent Inu 
different  Initial  temperature  and  pressure  conditions  in  the  driver  hi so. 

9  was  different  from  the  other  8  In  that  the  shock  tube  geometry  was  a 
little  different  (see  Table  1)  and  the  valve  was  opened  to  a  u:x;mum  of 
43.75%  instead  of  100%. 

The  valve  was  represented  as  a  parabolic  constriction  wHh  a  tots1 
length  of  2  meters,  located  In  the  middle  of  the  throat  The  valve 
opening  was  represented  as  a  linear  function  of  time,  from  0.1%  open  at 
time  zero  to  100%  open  (43.75%  open  for  case  9)  at  tre  specified  valve 
opening  time.  In  addition  to  these  calculations  w  in  tu  va!v.  ,;:e-  nj 
qradually,  a  reference  calculation  was  done  with  the  valve  opened 
instantaneously  for  each  case,  to  insure  adequate  resolution  of  thr  jho  -■ 
wave  structure  with  the  mesh  used  and,  in  some  instance'  to  compare 
certain  features  of  the  wave  structures. 
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Table  1:  Shock  Tube  Geometry 


»  m 


driver  length 
driver  diameter 
converging  nozzle  length 
throat  length 
throat  diameter 
diverging  nozzle  length 
expansion  tube  length 
expansion  tube  diameter 

(all  dimensions 


cases  1*6 

cases  7-8 

case  9 

80. 

20. 

20. 

5.490 

5.490 

5.440 

2. 

2. 

2. 

4. 

4. 

4. 

2.743 

2.743 

i.8: 

4. 

4. 

4. 

110. 

no. 

110. 

14.494 

14.494 

14.494 

are  in  meters) 
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TABLE  2:  Calculation  Results 


case 

pressure 

ratio 

temperature 

ratio 

nominal 

over¬ 

pressure 

[kPa] 

reference 

over¬ 

pressure 

[kPa] 

maximum 

valve 

opening 

time 

[second] 

1 

118.0 

2.235 

241.3 

247.0 

0.063 

2 

102.5 

2.009 

206.8 

213.6 

0.062 

3 

83.5 

1.825 

172.4 

174.5 

0.061 

4 

69.2 

1.652 

137.9 

142.7 

0.055 

5 

54.0 

1.471 

103.4 

107.2 

0.044 

6 

39.5 

1.228 

68.9 

74.0 

0.050 

7 

19.5 

1.094 

34.5 

39.6 

0.035 

8 

7.0 

1.040 

13.8 

15.57 

0.021 

9 

8.9 

1.040 

13.8 

4.634 

0.026 

* 
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The  calculations  were  done  on  INEL's  CRAY  X/MP-24  computer.  Two  minor 


modifications  were  made  to  the  BRL-Q1D  code:  the  sizes  of  the 
grid-dependent  arrays  were  Increased,  and  the  plot  and  print  output  of 
time-dependent  arrays  was  not  started  until  the  shock  reached  the  first 
measurement  station.  For  cases  B  and  9,  the  code  was  further  modified  by 
overlaying  some  large  arrays,  to  decrease  computer  memory  requirements  and 
thereby  Improve  computer  access  and  turn-around  time. 

The  BRL-Q10  code  Includes  two  options  for  the  numerical  method:  the 
Beam-Warming  and  the  MacCormack  methods.  For  these  calculations,  the 
Beam-Warming  method  was  used,  With  the  BRL- recommended  values  of  1.0  for 
the  Courant  number  and  0.04  for  the  smoothing  coefficient.  Preliminary 
estimates  suggested  that  10000  grid  points  would  be  adequate  for  these 
calculations,  and  this  number  was  used  successfully  for  cases  1-6.  More 
grid  points  were  required  for  cases  7-9. 
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3.  RESULTS 


The  results  for  case  1  are  given  here  In  some  detail;  results  for  the 
other  cases  are  given  In  much  less  detail,  only  as  Is  necessary  to 
describe  significant  differences  from  the  results  for  case  1. 

Figures  2-15  Indicate  the  temporal  and  spatial  dependence  of  some 
important  variables  calculated  by  the  BRL-Q1C  code  for  case  1  with  the 
valve  Instantaneously  opened  100%.  In  Figures  2-13,  the  horizontal  axis 
is  the  normalized  distance  along  the  length  of  the  shock  tube.  The  actual 
total  length  corresponding  to  the  normalized  length  of  1  is  indicated  by 
the  value  of  Lref  (200  meters  for  case  1)  printed  at  the  top  left  corner 
of  the  page.  The  profile  of  the  shock  tube  Is  Indicated  above  the  data 
o'ot,  on  the  same  horizontal  scale  as  the  data  plot.  The  dependent 
variables  plotted  on  the  vertical  scale  In  Figures  2-13  are  also 
normalized  to  reference  values  as  described  In  the  BRL-Q10  code 
documentation. The  number  printed  just  above  the  right  end  of  each 
curve  in  Figures  2-13  is  the  time  corresponding  to  that  curve,  In 
seconds.  In  Figures  10-13,  the  horizontal  scale  is  expanded  enough  so 
that  the  individual  points  of  the  spatial  grid  are  discernible,  Indicated 
by  the  vertical  lines  in  the  shock  tube  profile  drawn  above  the  graph 
proper.  These  graphs  with  high  spatial  resolution  clearly  indicate  the 
:patial  spread  of  the  computed  shock  front,  which  is  useful  in  evaluating 
the  adequacy  of  the  numerical  method. 


I 
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Figures  14  and  15  are  graphs  of  overpressure  (solid  line)  and  dynamic 
pressure  (dashed  line)  versus  time,  at  the  test  location  specified  by  the 
value  listed  for  X-sta  at  the  top  of  the  page.  The  Indicated  distance  of 
97.5  meters  from  the  beginning  of  the  expansion  tube  corresponds  to  the 
required  103.5  meters  from  the  center  of  the  valve. 

The  data  of  most  Immediate  Interest  Is  the  overpressure  versus  time 
data  Indicated  by  the  solid  curves  In  Figures  14  and  15.  This 
overpressure  rises  abruptly  --  In  a  fraction  of  a  millisecond  --  when  the 
shock  arrives  at  the  test  location.  At  the  top  of  this  step  Is  a  small 
overshoot  that  is  an  artifact  of  the  numerical  procedures  and  Is  not 
indicative  of  a  real  phenomenon;  the  actual  overpressure  step  is  better 
indicated  by  extrapolating  the  almost-horlzontal  curve  from  the  right  o^ 
the  step  into  the  region  of  the  step  itself.  This  extrapolation  procedure 
was  used  to  estimate  the  overpressure  step  In  all  the  studies  reported 
here. 

After  the  initial  shock  transient,  the  overpressure  approaches  a 
constant  value,  a  plateau  value,  which  Is  different  for  each  of  the  9 
cases.  This  plateau  value  is  used  as  the  reference  overpressure  for  each 
case. 
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Figures  16-29  Illustrate  iha  affect  of  the  valve  opening  gradually 
Initaad  of  Instantaneously.  These  figures  riprtsant  a  valve  opening  time 
of  65  milliseconds,  which  Is  close  to  the  maximum  allowable  valve  opening 
time  for  case  1.  The  major  effect  Is  to  decrease  the  magnitude  of  the 
overpressure  step  and  to  Increase  the  slope  of  the  overpressure  versus 
time  curve  after  the  Initial  shock.  The  rise  time  for  the  Initial  step  Is 
not  noticeably  changed  by  the  gradual  opening  of  the  valve.  Since  the 
slope  of  the  overpressure  versus  time  curve  after  the  Initial  shock  Is 
practically  zero  on  a  millisecond  time  scale,  the  acceptance  criterion  Is 
essentially  equivalent  to  the  requirement  that  the  magnitude  of  the 
overpressure  step  with  the  valve  opened  gradually  be  at  least  95%  of  the 
reference  overpressure.  Thus,  the  procedure  for  estimating  the  maximum 
allowable  valve  opening  time  was  to  try  different  opening  times  and  check 
whether  the  shock  overpressure  step  was  at  least  95%  of  the  reference 
value.  (Of  course  the  shock  rise  time  was  also  checked  for  each  run,  to 
be  sure  that  It  was  less  than  0.5  ms.) 

This  case  1  description  Is  representative  of  the  results  for  cases 
1-5.  Case  6  and  7  gave  results  similar  to  those  described  for  case  1, 
except  that  the  overpressure  decreased  slightly  Instead  of  Increasing 
after  the  Initial  shock  for  the  calculations  In  which  the  valve  was  opened 
instantaneously.  This  difference  was  not  Important  in  estimating  the 
maximum  acceptable  valve  opening  time. 


Cases  8  and  9  gave  qualitatively  different  results,  as  Indicated  In 
Figures  30-39  for  case  8  with  the  valve  Instantaneously  opened.  In  these 
cases,  the  Initial  shock  was  followed  after  a  few  milliseconds  by  a  second 
shock  that  raised  the  overpressure  to  a  value  substantially  higher  than 
that  of  the  first  shock.  The  timing  and  magnitude  of  the  second  shock 
depend  on  the  details  of  the  geometry  In  the  throat  region  of  the  shock 
tube.  Figures  40-93  show  some  detailed  plots  for  case  8  with  the  valve 
Instantaneously  opened,  and  Figure  94  shows  an  x-t  plot  for  this  case. 

The  second  shock  moves  a  little  faster  that  the  first  shock,  but  does  not 
quite  overtake  the  first  shock  In  the  103.5  meter  travel  distance  to  the 
test  station.  It  appeared  that  moving  the  valve  closer  to  the  end  of  the 
throat  would  decrease  the  initial  delay  between  the  first  and  second 
shocks  and  perhaps  allow  the  second  shock  to  overtake  and  merge  with  the 
first  shock  before  they  reached  the  test  location.  This  change  In  the 
valve  location  was  tried  and  found  effective  for  this  particular  problem; 
unfortunately,  this  created  another  shock  even  later,  so  that  the  overall 
effect  was  not  significantly  changed. 

The  general  conclusion  drawn  from  this  examination  of  case  8  Is  that 
the  details  of  the  shock  tube  geometry  near  the  valve  can  cause  Irregular 
behavior  of  the  overpressure  versus  time  plot  for  the  shock  front,  more 
complicated  than  the  simple  step  followed  by  an  almost  flat  region.  In 
fact,  this  same  irregular  behavior  occurs  in  all  cases  (note  Figure  3), 
but  In  the  higher  pressure  cases  the  Irregular  structure  Is  merged  Into 
one  simple  shock  before  the  shock  reaches  the  measurement  station. 
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Of  course  the  simple  shock  tube  geometry  assumed  here  arid  the 
one-dimensional  fluid  flow  calculation  used  here  may  not  accurately 
represent  the  actual  shock  tube  valve,  throat,  and  nozzle  behavior; 
nevertheless,  It  would  not  be  surprising  If  phenomena  qualitatively 
similar  to  the  double  shock  behavior  calculated  for  case  8  actually 
occurred  In  the  real  L8/TS.  The  question  of  Immediate  Importance  here  Is 
whether  the  model  Is  realistic  or  Informative  enough  to  be  useful,  and 
whether  (and  how)  to  proceed  with  the  studies  of  acceptable  valve  opening 
times  for  cases  8  and  9.  The  BRL  personnel  consulted^3)  on  this  matter 
suggested  prodeedlng  with  the  calculations,  with  the  understanding  that 
the  results  might  not  be  taken  at  face  value.  Thus,  the  results  presented 
here  for  cases  8  and  9  should  be  regarded  as  plausible  but  perhaps  not 
very  accurate  estimates  of  the  actual  shock  tube  behavior.  These  results 
are  a  little  sketchy,  but  more  detailed  studies  did  not  seem  warranted  In 
view  of  the  high  cost  and  limited  validity  of  computer  calculations  for 
the  low-pressure  cases  like  cases  8  and  9. 

Overpressure  versus  time  plots  for  cases  8  and  9  are  shown  In  Figure 
38  and  Figures  95-100.  The  overpressure  step  associated  with  the  first 
shock  front  appears  to  vary  In  a  consistent  manner  with  valve  opening 
time,  if  the  valve  Is  opened  gradually  Instead  of  Instantaneously,  and  the 
results  obtained  from  considering  only  this  first  shock  front  fit 
reasonably  well  with  the  results  from  cases  1-7.  Therefore,  despite  the 
irregular  nature  of  the  overpressure  versus  time  curve  near  the  shock 
front,  cases  8  and  9  were  finally  treated  In  the  same  manner  as  the  other 
cases,  ..with  the  overpressure  step  at  the  first  shock  front  being 
considered  the  critical  factor. 
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For  the  case  9  calculations  with  valve  opening  tines  of  20  and  30 
milliseconds,  the  shock  front  rise  time  was  a  little  greater  than  the 
required  0.5  milliseconds.  However,  previous  experience  indicates  that 
the  parameters  of  Interest  are  not  very  sensitive  to  this  factor,  and  the 
high  calculation  cost  prohibited  repeating  these  calculations  with  smaller 
step  sizes  to  obtain  the  required  rise  time. 

The  results  of  these  valve  opening  time  studies  are  summarized  In 
Figures  101-103,  which  Indicate  the  magnitude  of  the  shock  overpressure 
step  versus  valve  opening  time  for  each  of  the  9  cases  studied. 
Interpolation  of  these  curves  gives  the  maximum  allowable  valve  opening 
time  for  each  case,  listed  in  Table  2  and  graphed  versus  driver  pressure 
to  ambient  pressure  ratio  In  Figure  104.  This  curve  does  show  a  genera! 
trend,  but  It  Is  not  smooth.  This  lack  of  smoothness  is  at  least  partly 
attributed  to  the  fact  that  the  maximum  allowable  valve  opening  time  Is 
not  a  function  of  pressure  ratio  alone;  It  also  depends  on  the  temperature 
ratio,  which  is  not  a  perfectly  smooth  function  of  the  pressure  ratio  (see 
Figure  105)  for  the  cases  considered  here. 


4.  CONCLUSIONS 


Th«  calculations  don*  here  Indicate  that  tha  maximum  allowable  valve 
opening  time  It  about  63  mlllltecondi  for  cate  1,  and  tomewhat  lett  for 
the  other,  lower-pressure  catet,  at  Indicated  In  Table  2.  The  retults  for 
catet  1-7  are  believed  to  be  real  1 ttlc;  the  retults  for  cases  8  and  9  must 
be  regarded  with  some  skepticism.  Mora  reliable  results  would  require  a 
multidimensional  computer  model  that  Is  somewhat  more  complicated  and 
substantially  more  expensive  to  run  than  the  model  used  In  this  work. 
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Figure  1.  Oi'agram  (not  to  scale)  of  shock  tube  for  Case 
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Figure  3.  Magnified  pressure  versus  distance,  case  1, 
zero  opening  time 


Figure  4.  Density  versus  distance,  case  1,  zero  opening  time 
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Figure  6.  Velocity  versus  distance,  case  1,  zero  opening  time 
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Figure  7.  Magnified  velocity  versus  distance,  case  1, 
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Figure  10.  Shock  front  detail,  pressure  versus  distance, 
case  1,  zero  opening  time 
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Figure  11.  Shock  front  detail,  density  versus  distance, 
case  1,  zero  opening  time 
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Figure  12.  Shock  front  detail,  velocity  versus  distance, 
case  1,  zero  opening  time 


PLOT  12  08.19.40  THUR  10  DEC,  1987  J08-PI03QD86,  ISSCO  OlSSPLfl  10.0 


Figure  13.  Shock  front  detail,  sound  speed  versus  distance, 
case  1,  zero  opening  time 
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Figure  18.  Density  versus  distance,  case  1,  65ms  opening  time 
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Figure  19.  Magnified  density  versus  distance,  case  1, 
65ms  opening  time 
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Figure  25.  Shock  front  detail,  density  versus  distance, 
case  1,  65ms  opening  time 
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Figure  26. 


Shock  front  detail,  velocity  versus  distance, 
case  1,  65ms  opening  time 
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Figure  27.  Shock  front  detail,  sound  speed  versus  distance, 
case  1,  65ms  opening  time 


STATIC  &  DYNAMIC(Z>cts7i)  PRESSURE,  kPa 

25  50  75  100  125  150  175  200  225  £ 


SHOCK  TUBE  WITH  CONICAL  DRIVER  AND  NOZZLE 


CASE: 

=  lu:  BRL1 

P-so 

=  24-1.2  kPa 

(  34. ■>(  psi) 

L-ref 

=  300.0  m 

X-sta  =  97.50 

m 

t-a 

=  204.8  ms 

L-drv 

=  80.00  m 

P  -drv  =  118.0 

atm 

PPD 

=  0.015  s 

(  o.trvi  s) 

V-drv 

=  1929.  m3 

P  a  mb  =  101.3 

kPa 

I  -so 

=  3.502  kPa-s 

(  0.01'  k'i) 

L-dvn 

=  110.0  m 

T-arnb  =  288.2 

K 

Q-s 

=  153.2  kPa 

L-rwc 

=  0.000  m 

T4/T1  =2.235 

I-dy  n 

=  2.224  kPa-'S 

i  o  or-:  k  r  i 

PRESSURE  TIME  HISTORY 


TIME,  3 


Figure  28.  Pressure  versus  time,  case  1,  65ms  opening  time 
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Figure  29.  Detail,  pressure  versus  time,  case  1,  65ms  opening  time 
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Figure  31.  Magnified  pressure  versus  distance,  case  8, 
zero  opening  time 
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Figure  35.  Magnified  velocity  versus  distance,  case  8 
zero  opening  time 
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Figure  39.  Detail,  pressure  versus  time,  case  8,  zero  opening  time 
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Figure  40.  Detail  study,  pressure  versus  distance,  part  1 
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Figure  47.  Detail  study,  density  versus  distance,  part  2 
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Figure  48.  Oetail  study,  density  versus  distance,  part  3 
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Figure  49.  Detail  study,  magnified  density  versus  distance,  part  1 
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Figure  54.  Detail  study,  velocity  versus  distance,  part  3 
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Figure  55.  Oetall  study,  magnified  velocity  versus  distance,  part  1 
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Figure  61.  Detail  study,  magnified  sound  speed  versus  distance, 
part  1 
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Figure  63.  Detail  study,  magnified  sound  speed  versus  distance, 
part  3 
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Figure  65.  Detail  study,  pressure  versus  distance 
0.3  to  0.7,  part  2 
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Figure  66.  Detail  study,  magnified  pressure  versus  distance 
0.3  to  0.7,  part  1 
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Figure  67.  Detail  study,  magnified  pressure  versus  distance 
0.3  to  0.7,  part  2 
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Figure  70.  Detail  study,  magnified  density  versus  distance 
0.3  to  0.7,  part  1 
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Figure  71.  Detail  study,  magnified  density  versus  distance 
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Figure  73.  Detail  study,  velocity  versus  distance 
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Figure  80.  Detail  study,  pressure  versus  distance  0.38  to  0.48 
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Figure  81.  Detail  study,  density  versus  distance  0.38  to  0.48 
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Figure  82.  Oetail  study,  velocity  versus  distance  0.38  to  0.48 
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Figure  83.  Detail  study,  sound  speed  versus  distance  0.38  to  0.48 
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Figure  84.  Detail  study,  pressure  versus  distance  0.39  to  0.43 
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Figure  85.  Detail  study,  density  versus  distance  0.39  to  0.43 
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Figure  87.  Detail  study,  sound  speed  versus  distance  0.39  to  0.43 
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Figure  88.  Detail  study,  pressure  versus  time  at  -4  meters 
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Figure  89.  Detail  study,  pressure  versus  time  at  -2 
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Figure  92.  Detail  study,  pressure  versus  time  at  4  meters 
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Figure  95.  Pressure  versus  time,  case  8,  20ms  orsr.iiq 
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TEMPERATURE  RATIO 


Figure  105.  Driver  temperature  ratio  versus  driver  pressure  ratio 
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